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' Abstract 

o '. 

\^ . We review the status of the QCD corrected cross sections and kinematic 

O I distributions for the production of a Higgs boson in association with top quark 

• or bottom quark pairs at the Fermilab Tevatron and at the LHC. Results for 

bbH production are presented in the Minimal Supersymmetric Model, where 
Ph! the rates can be greatly enhanced relative to the Standard Model rates. We 

place particular emphasis on theoretical uncertainties due to renormalization 
and factorization scale dependence and on the uncertainties coming from the 
Parton Distribution Functions. 



1 Introduction 

A light Higgs boson is preferred by precision fits of the Standard Model (SM) and 
also theoretically required by the Minimal Supersymmetric extension of the Standard 
Model (MSSM). The production of a Higgs boson in association with a heavy quark 
and antiquark pair, both ti and bb, at the Tevatron and the Large Hadron Collider 
(LHC) will be sensitive to the Higgs-fermion couplings and can help discriminate 
between models. 
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2 RESULTS FOR TTH PRODUCTION 



The associated production of a Higgs boson with a pair of tt quarks has a distinc- 
tive signature and can give a direct measurement of the top quark Yukawa couphng. 
This process is probably not observable at the Tevatron, but will be a discovery chan- 
nel at the LHC for M/j < 130 GeV. The associated production of a Higgs boson with 
a pair of bb quarks has a small cross section in the Standard Model, and can be used 
to test the hypothesis of enhanced bottom quark Yukawa couphngs in the MSSM 
with large values of tan/?. Both the Tevatron and the LHC will be able to search 
for enhanced bbh production, looking for a final state containing no bottom quarks 
(inclusive production), one bottom quark (semi-inclusive production) or two bottom 
quarks (exclusive production). 

The rates for tih production at the Tevatron and the LHC have been calculated 
at NLO QCD several years ago[l, 2, 3, 4, 5, 6]. The theoretical predictions for bbh 
production at hadron colliders involve several subtle issues, and depend on the number 
of bottom quarks identified in the final state. In the case of no or only one tagged 
bottom quark there are two approaches available for calculating the cross sections for 
bbh production, called the four flavor number schemes (4FNS)[7, 8] and flve flavor 
number scheme (5FNS)[9]. The main difference between these two approaches is that 
the 4FNS is a fixed-order calculation of QCD corrections to the gg and gg-induced bbh 
production processes, while in the 5FNS the leading processes arise from bg {bg) and 
U) initial states and large coUinear logarithms are resummed using a pertubatively 
defined bottom quark Parton Distribution Function (PDF). Very good agreement is 
found for the NLO QCD corrected cross sections for bb Higgs associated production 
when the two schemes are compared[10, 11]. 

In the following sections, we present numerical results at NLO QCD for tih and 
bbh production at the Tevatron and the LHC. If not stated otherwise, numerical 
results have been obtained in the 4FNS. We emphasize theoretical uncertainties from 
scale and PDF uncertainties and also present differential cross sections at NLO for 
Wh production in the case when two b quarks are tagged. 

2 Results for tih Production 

The observation of a tih final state will allow for the measurement of the tih Yukawa 
couphng. If Mh < 130 GeV, pp — > tih is an important discovery channel for a SM- 
hke Higgs boson at the LHC (-\/s = 14 TeV) [12, 13]. Given the statistics expected 
at the LHC, pp tth, with h bb,T~^T~ ,W~^W~ will be instrumental for the 
determination of the couplings of the Higgs boson. Precisions of the order of 10- 
15% on the measurement of the top quark Yukawa couphng can be obtained with 
integrated luminosities of 100 fb~^ per detector[14, 15, 16, 17]. 

The impact of NLO QCD corrections on the total cross section for pp, pp — > tih 
production in the Standard Model is illustrated in Fig. 1[6, 5, 1, 3] and Fig. 2[6, 3]. 
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3 RESULTS FOR BBH PRODUCTION 
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Figure 1: Total LO and NLO cross sections for pp,pp 
v/i= 1.96 TeV and ^/s = 14 TeV, for fi=mt + Mh/2. 



tth as functions of Mh, at 



The dependence of the total cross sections on the renormalization and factorization 
scales is strongly reduced at NLO as shown in Fig. 2. The numerical results at 
NLO are obtained using CTEQ4M (Fig. 2 (l.h.s.)), CTEQ5M (Fig. 2 (r.h.s.)), and 
CTEQGIvI (Fig. 1) parton distribution functions. The NLO cross section is evaluated 
using the 2-loop evolution of as(/i) with af^'^{Mz) = 0.116 (Fig. 2 (l.h.s.)) and 
a^^^{Mz) = 0.118 (Fig. 2 (r.h.s.)) and Fig. 1), and mt = 174 GeV. The renormal- 
ization/factorization scale dependence, uncertainty on the PDFs, and the error on 
the top quark pole mass, mt, are estimated to give a 15-20% uncertainty. 



3 Results for bbh Production 

The bbh production processes are only relevant discovery modes in the IvlSSIvl with 
large tan (3. To a good approximation, the predictions for the IVISSIVI rates can easily 
be derived from the Standard Ivlodel results by rescaling the Yukawa couplings [10]. 
The dominant IMSSIVI radiative correction to bbh production can be taken into account 
by including the IvISSIvI corrections to the bbh vertex only, i.e. by replacing the tree 
level Yukawa couplings by the radiative corrected ones. We follow the treatment of 
the program FeynHiggs [18, 19] and take into account the leading, tan/3 enhanced, 
radiative corrections that are generated by gluino-sbottom and chargino-stop loops. 
For large tan (3, the bottom quark Yukawa coupling is enhanced and the top quark 
Yukawa coupling coupling is strongly suppressed, resulting in a IvlSSIvI bbh cross 
section that is about three orders of magnitude larger than the Standard IVlodel cross 
section. For the Tevatron, we calculate the production rates for the lightest IvISSIM 
Higgs boson, h^, while for the LHC we consider the rate for the heavier neutral Higgs 
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3.1 Total Cross Sections for bbh Produc13onRESULTS FOR BBH PRODUCTION 
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Figure 2: Dependence of (yLo,NLo{pPiPP ~^ tth) on the renormalization/factorization 
scale /i, at ^/s = 2 TeV (l.h.s.jand yi=14 TeV (r.h.s.), for Mh = l20 GeV. 



boson, H^} 

In the numerical evaluation of cross sections for the exclusive and semi-inclusive 
channels {bbh and bh + bh production), it is required that the final state bottom 
quarks have pt>20 GeV and pseudorapidity |?7|<2.0 for the Tevatron and |?7|<2.5 
for the LHC. In the NLO real gluon emission contributions, the final state gluon 
and bottom quarks are considered as separate particles only if their separation in the 
pseudorapidity-azimuthal angle plane, AR = a/ (Ar])^ + (A^)^, is larger than 0.4. For 
smaller values of AR, the four momentum vectors of the two particles are combined 
into an effective bottom/anti-bottom quark momentum four- vector. 

If not stated otherwise, the numerical results at NLO are obtained using CTEQ6M 
PDFs, the 2-loop evolution of ^^(/x) with a™(Mz) = 0.118, and the WS renormal- 
ization scheme for the bottom quark mass and Yukawa coupling with 2-loop renor- 
malization group improved MS masses. The bottom quark pole mass is chosen to be 
rub = 4.62 GeV. 

3.1 Total Cross Sections for bbh Production 

We present total cross section results at NLO in the 4FNS in Fig. 3 for associated bb 
Higgs production in the MSSM with tan /3 = 40. The bands represent the theoretical 
uncertainty due to the residual scale dependence. They have been obtained by varying 

^We assume Msusy = 1 TeV, Mg = 1 TeV, At ^ At ^ 2 TeV Ab = At = 25 GeV {H°), 
^ = M2 = 200 GeV and ^ = Afz = 1 TeV (H^). For M^,o = 120 GeV, the bbh° coupling is 

enhanced by a factor of 33 relative to the SM coupling, while for M^o between 200 and 800 GeV, 
the coupling is enhanced by a factor of 27 relative to the SM coupling. 
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3.1 Total Cross Sections for bbh Produc13onRESULTS FOR BBH PRODUCTION 
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Figure 3: Total cross sections for pp,pp — > bbh in the MSSM in the 4FNS at NLO for 
the Tevatron and the LHC in the MSSM with tan/5 = 40 and with 0,1 or 2 6 quarks 
identified. The Tevatron (LHC) plot is for the lightest (heaviest) neutral Higgs boson, 
(H^)- The error bands have been obtained by varying the renormalization and 
factorization scales as described in the text. 



the renormalization {fir) and factorization (fif) scales independently from /io/4 to fiQ, 
where fio = mb + Mh/2. 

If the outgoing bottom quarks cannot be observed then the dominant MSSM 
Higgs production process at large tan/3 is gg —>■ {bb)h (the curve labelled '0 b'). 
The inclusive cross section is experimentally relevant only if the Higgs boson can be 
detected above the background without tagging bottom quarks. At the LHC, this 
process can be identified at large tan (5 by the decays to and r+r^ for the 

heavy Higgs bosons, and , of the MSSM. At the Tevatron this process, with 
T^T~ , has been used to search for the neutral MSSM Higgs boson. If a single 
bottom quark is tagged then the final state is bh or bh (the curve labelled '1 b'). 
Although requiring a b quark in the final state significantly reduces the rate, it also 
reduces the background. A recent Tevatron study [20] used the search for neutral 
MSSM Higgs bosons in events with three bottom quarks in the final state {bh^ + bh^ 
production with — > bb) to impose limits on the tan (3 and M^o parameter space. 

Finally, we show the fully exclusive cross sections for bbh production, where both 
the outgoing b and b quarks are identified (the curve labelled '2 b'). The exclusive 
measurement corresponds to the smallest cross section, but it also has a significantly 
reduced background. Moreover, both the exclusive and semi-inclusive bbh production 
modes are the only ones that can unambiguously measure the bottom quark Yukawa 
coupling. 
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3.2 Differential Cross Sections for bbh PAdiRmNLTS FOR BBH PRODUCTION 
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Figure 4: The relative corrections daNLo/do^Lo — 1 for the pr distribution of the 
bottom or anti-bottom quark with the largest pt {pT,max) (left) and of the SM Higgs 
boson (p^) (right) to bbh production in the SM at the Tevatron (with y/s = 2 TeV 
and /i = 2/io) and the LHC (with y^=14 TeV and /i = 4/io). 




3.2 Differential Cross Sections for bbh Production 

In assessing the impact of the NLO corrections it is particularly interesting to study 
the kinematic distributions. In Figs. 4 and 5 we illustrate the impact of NLO QCD 
corrections on the transverse momentum and pseudorapidity distribution of the SM 
Higgs boson and the bottom quark by showing the relative correction, daNLO I daio — ^ 
(in percent) for the exclusive case {hbh where both b quarks are observed). For the 
renormalization/factorization scale we choose /i = 2/io at the Tevatron and /z = 4/io at 
the LHC, with = mh + Mh/2, and use the CTEQ5 set of PDFs. As can be seen, the 
NLO QCD corrections can considerably affect the shape of kinematic distributions, 
and their effect cannot be obtained from simply rescaling the LO distributions with 
a K-factor of (JiVLo/<7Lo = 1-38 ± 0.02 (Tevatron, /i = 2/io) and cr^r^o/<7Lo = l-H i 0.03 
(LHC, /i = 4/io).' 

3.3 PDF and Renormalization/Factorization Scale Uncertain- 
ties 

A major source of theoretical uncertainty for cross section predictions comes from 
the PDFs. We study the uncertainties of semi-inclusive bh production rates from 
the uncertainties in the PDFs using the CTEQ PDF sets [21]. First, the central 
value cross section ao is calculated using the global minimum PDF (i.e. CTEQ6M). 
The calculation of the cross section is then performed with the additional 40 sets 
of PDFs to produce 40 different predictions, cxj. For each of these, the deviation 
from the central value is calculated to be Aa^^ = \ai — (TqI when cTi^ao- Finally, to 

^The kinematic distributions have been calculated within the Standard Model and using the on- 
shell scheme for the definition of the b quark mass, but we see a similar behavior when using the 
MS bottom quark Yukawa coupling. 
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3.3 PDF and Renormalization/FactoriziAioMESMJU&cESlMiMBsH PRODUCTION 




Figure 5: The relative corrections daMLo/d(^Lo ~ 1 for the t] distribution of the bottom 
quark 77^ (left) and of the SM Higgs boson (jih) (right) to bbh production in the SM 
at the Tevatron (with ^/s = 2 TeV and /i = 2yUo) and the LHC (with yi=14 TeV and 
/i = 4/io). 
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Figure 6: Normalized cross sections for Higgs production with one h jet at the Teva- 
tron (l.h.s) and the LHC (r.h.s) showing the uncertainty from PDFs for both the gg 
(4FNS) and hg (5FNS) initial states. 



obtain the uncertainties due to the PDFs the deviations are summed quadratically 
as Acr^ = \l^ji^^f^ and the cross section including the theoretical uncertainties 

arising from the PDFs is quoted as ctqI^^^-. 

In Fig. 6, we plot the normalized total SM NLO cross sections for semi-inclusive 
hh production, calculated in the 5FNS {hg —>■ bh) as implemented in MCFM [22] and 
in the 4FNS [gg b{b)h), and compare their respective uncertainties due to the 
PDFs. We see that, at both the Tevatron and the LHC, the PDF uncertainties are 
almost identical for both the gg and bg initial states. 

In Figs. 7 and 8 we compare the uncertainties from residual scale dependence and 
the PDFs on the example oibg bh (5FNS) at the Tevatron and LHC respectively [10]. 
Here, we perform the comparison for both the total cross section (left) and the to- 
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4 CONCLUSION 




MJGeVl Mj(GeV) 

Figure 7: Comparison between theoretical uncertainties due to scale dependence and 
uncertainties arising from the PDFs at the Tevatron for semi-inclusive bh production 
in the Standard Model. In the right-hand plot, both uncertainty bands have been 
normalized to the central value of the total cross section ctq. 



tal cross section normalized to the central value calculated with CTEQ6M (right). 
Similar results are obtained in the 4FNS. 

From Fig. 8 one can see that, at the LHC, the theoretical uncertainty is dominated 
by the residual scale dependence. Due to the large center of mass (c.o.m.) energy of 
the LHC, the gluons and bottom quarks in the initial state have small momentum 
fraction (x) values and, hence, small PDF uncertainties typically in the 5-10% range. 

In contrast, due to the smaller c.o.m. energy, the PDF uncertainties at the Teva- 
tron (Fig. 7) are comparable and even larger than the uncertainties due to residual 
scale dependence over the full Higgs mass range. The smaller c.o.m. energy results 
in higher-x gluons and bottom quarks in the initial state which corresponds to large 
PDF uncertainties in the 10-30% range. 

4 Conclusion 

The NLO cross sections for tih and bbh have been presented for the Tevatron and the 
LHC with emphasis on the renormalization/factorization scale and PDF dependences. 
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1 Hadronic Higgs Production with Heavy Quarks 
at the Tevatron and the LHC 



1.1 Introduction 

A light Higgs boson is preferred by precision fits of the Standard Model (SM) and 
also theoretically required by the Minimal Supersymmetric extension of the Standard 
Model (MSSM). Searches at both the Tevatron and the Large Hadron Collider (LHC) 
will play a crucial role in discriminating between models. In this context, the pro- 
duction of a Higgs boson in association with a heavy quark and antiquark pair, both 
tt and bb, plays a very important role. 

The associated production of a Higgs boson with a pair of ti quarks has a distinc- 
tive signature, and can give a direct measurement of the top quark Yukawa coupling. 
This process is probably not observable at the Tevatron, but will be a discovery chan- 
nel at the LHC for < 130 GeV. The associated production of a Higgs boson with 
a pair of bb quarks has a small cross section in the Standard Model, and can be used 
to test the hypothesis of enhanced bottom quark Yukawa couplings in the MSSM 
with large values of tan/?. Both the Tevatron and the LHC will be able to search 
for enhanced bbh production, looking for a final state containing no bottom quarks 
(inclusive production), one bottom quark (semi-inclusive production) or two bottom 
quarks (exclusive production). 

The theoretical prediction for the bbh production rate at hadron colliders involves 
several subtle issues, and depends on the number of bottom quarks identified in the 
final state. In the case of no or only one tagged bottom quark there arc two approaches 
available for calculating the cross sections for bbh production, called the four fiavor 
(4FNS) and five flavor (5FNS) number schemes. The main difference between these 
two approaches is that the 4FNS is a fixed-order calculation of QCD corrections to the 
gg and gg-induced bbh production processes, while in the 5FNS the leading processes 
arise from bg (bg) and bb initial states and large collinear logarithms are resummed 
using a pertubatively defined bottom quark Parton Distribution Function (PDF). The 
NLO QCD corrected cross sections for b Higgs associated production are in excellent 
agreement in when the two schemes are compared. 

In the following sections, we present NLO numerical results for tih and bbh pro- 
duction at the Tevatron and the LHC. We emphasize theoretical uncertainties from 
scale and PDF uncertainties and also present differential cross sections at NLO. 

1.2 Results for tih Production 

The observation of a tih final state will allow for the measurement of the tih Yukawa 
coupling. Observing pp — > tih at the Tevatron (y^= 1.96 TeV) will require very high 
luminosity and will probably be beyond the machine capabilities. On the other hand. 
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Figure 1: <Jnlo{pp tih) and a^oipv ~^ tth) as functions of Mh, at 1/5 = 14 TeV, for 
fi = mt + Mh/2 and /x = 2mj + Mh- 

if Mh < 130 GeV, pp tth is an important discovery channel for a SM-like Higgs 
boson at the LHC (-^5 = 14 TeV) [???]. Given the statistics expected at the LHC, 
pp — > tih, with h bb, t^t~ , W~^W~, '-ff will be instrumental for the determination 
of the couplings of the Higgs boson [?????]. Several studies show that precisions 
of the order of 10-15% on the measurement of the top quark Yukawa coupling can be 
obtained with integrated luminosities of 100 fb~^ per detector. 

The impact of NLO QCD corrections on the total cross section for pp — > tth 
production in the Standard Model is illustrated in Fig. l.[? ? ? ? ] The numerical 
results are obtained using CTEQ6M parton distribution functions. The NLO cross 
section is evaluated using the 2-loop evolution of as{^) with a^^'-^^Mz) = 0.118 and 
Mt = 172 GeV. The renormalization/factorization scale dependence is estimated to 
give a ?-?% uncertainty. 

1.3 bbh Production 

The bbh processes are only relevant discovery modes in the MSSM with large tan f3. 
The predictions for the MSSM rates can easily be derived from the Standard Model 
results by rescaling the Yukawa couplings. The dominant MSSM radiative correction 
to bbh production can be taken into account by including the MSSM corrections to 
the bbh vertex only, i.e. by replacing the tree level Yukawa couplings by the radiative 
corrected ones. We follow the treatment of the program FeynHiggs [? ] and take 
into account the leading, tan (3 enhanced, radiative corrections that are generated by 
gluino-sbottom and chargino-stop loops. For large tan f3, the bottom Yukawa coupling 
is enhanced and the top Yukawa coupling coupling is strongly suppressed, resulting 
in a MSSM bbh cross section that is about three orders of magnitude larger than the 
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Standard Model cross section. For the Tevatron, we calculate the production rates 
for the lightest MSSM Higgs boson, h^, while for the LHC we consider the rate for 
the heavier neutral Higgs boson, H^} 

Numerical Results for bbh 

For the exclusive and semi- inclusive channels {bbh and bh + bh production), it is 
required that the final state bottom quarks have pr > 20 GeV and pseudorapidity 
1 77 1 < 2.0 for the Tevatron and |r7|<2.5 for the LHC. In the NLO real gluon emission 
contributions, the final state gluon and bottom quarks are considered as separate 
particles only if their separation in the pseudorapidity- azimuthal angle plane, AR = 



2 -|- (A0)2, is larger than 0.4. For smaller values of AR, the four momentum 



vectors of the two particles are combined into an effective bottom/anti-bottom quark 
momentum four- vector. 

Total cross sections for bbh production 

We present total cross section results at NLO in the 4FNS in Fig. 2 for associated 
b quark Higgs production in the MSSM with tan/? = 40. The bands represent the 
theoretical uncertainty due to the residual scale dependence. They have been obtained 
by varying the renormalization (/x,.) and factorization {fj,f ) scales independently from 
to /Iq, where iio — nib + Mfi/2. 
If the outgoing bottom quarks cannot be observed then the dominant MSSM 
Higgs production process at large tan/? is gg — > {bb)h (the curve labelled '0 b'). The 
exclusive measurement corresponds to the smallest cross section, but it also has a 
significantly reduced background. At the LHC, this process can be identified by the 
decays to and t+t~ for the heavy Higgs bosons, and of the MSSM. At 

the Tevatron this process, with /i° — > t+t~, has been used to search for the neutral 
MSSM Higgs boson. 

If a single bottom quark is tagged then the final state is bh or bh (the curve 
labelled 'lb'). A recent Tevatron study [? ] used the search for neutral MSSM Higgs 
bosons in events with three bottom quarks in the final state -|- production 
with /i° — > bb) to impose limits on the tan (3 and M40 parameter space. 

Finally, we show the fully exclusive cross sections for bbh production, where both 
the outgoing b and b quarks are identified (the curve labelled '2 b'). 



^We assume Msusy = 1 TeV, Mg = 1 TeV, At ^ At ^ 2 TeV {h"), Ab ^ At = 25 GeV (H"), 
^ = M2 = 200 GeV (/i°), and At = M2 = 1 TeV For M^o = 120 GeV, the bbh^ coupling is 

enhanced by a factor of 33 relative to the SM couphng, while for Mfjo between 200 and 800 GeV, 
the bbH^ coupling is enhanced by a factor of 27 relative to the SM coupling. 
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Figure 2: Total cross sections for pp,pp — > {bb)h in the MSSM in the 4FNS at NLO 
for the Tevatron and the LHC in the MSSM with tan f3 = 40. The Tevatron (LHC) 
plot is for the lightest (heaviest) neutral Higgs boson, (H^)- The error bands have 
been obtained by varying the renormalization and factorization scales as described in 
the text. 



Differential Cross Sections for Associated b Higgs Production 

In comparing the four and five flavor number schemes it is particularly interesting 
to compare the kinematic distributions. In Figs. 3-4, we compare the results for the 
NLO transverse momentum and pseudorapidity distributions of the light and heavy 
MSSM Higgs bosons, and H^, in both the 4FNS and 5FNS, at the Tevatron and 
the LHC. We see, in general, good agreement between the two schemes within their 
respective theoretical uncertainties, except in regions of kinematic boundaries. This is 
particularly dramatic in thep^ distributions in the 5FNS where, around j9^~ 20 GeV, 
a kinematic threshold induced in bg — > bh by the cut on the pt of the bottom quark 
causes the 5FNS NLO calculation to be unreliable. 

Finally, in Figs. 5 and 6 we illustrate the impact of NLO QCD corrections on the 
transverse momentum and pseudorapidity distribution of the SM Higgs boson and 
the bottom quark by showing the relative correction, daMiol da^o ~ 1 (in percent) 
for the exclusive case {hbh where both b quarks are observed). For the renormal- 
ization/factorization scale we choose ^ = 2fiQ at the Tevatron and n = at the 
LHC. As can be seen, the NLO QCD corrections can considerably affect the shape 
of kinematic distributions, and their effect cannot be obtained from simply rescaling 
the LO distributions with a K-factor of a f^i^o jo^o = 1-38 ± 0.02 (Tevatron, /i = 2/io) 
and (Ti^Lo/ ^Lo = 1-11 i 0.03 (LHC, /i = 4/io). The kinematic distributions have been 
calculated within the Standard Model and using the on-shell scheme for the definition 
of the b quark mass, but we see a similar behavior when using the MS bottom quark 
Yukawa coupling. 
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Figure 3: da/dp!^ in the MSSM at the Tevatron and the LHC for M;,o,jyo = 120 GeV 
and yUf. = yUi? = /io/2 for single 6-tag events. We show the NLO results in the 4FNS 
(solid) and 5FNS (dashed), using two different bin sizes, 2 GeV (left) and 12 GeV 
(right). 
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Figure 4: da/dr]h in the MSSM at the Tevatron and the LHC for M^o^ho = 120 GeV 
and /ir = /i/ = /^o/2 for single 6-tag events. We show the NLO results in the 4FNS 
(solid) and 5FNS (dashed). 
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Figure 5: The relative corrections daf^i^o/dcLo — 1 foi' the px distribution of the 
bottom or anti-bottom quark with the largest px {PT,max) (left) and of the SM Higgs 
boson (pt^) (right) to bbh production in the SM at the Tevatron (with ^/s = 2 TeV 
and /i = 2/io) and the LHC (with ^^=14 TeV and ^ = A^q). 




Figure 6: The relative corrections dapf^o/d'^Lo ~ 1 for the 1] distribution of the bottom 
quark rji, (left) and of the SM Higgs boson {rjh) (right) to bbh production in the SM 
at the Tevatron (with ^/s = 2 TeV and /i = 2/io) and the LHC (with ^/s = 14: TeV and 
/i = 4/io). 



PDF and Renormalization/Factorization Scale Uncertainties 

A major source of theoretical uncertainty for cross section predictions comes from the 
PDFs. We study the uncertainties of semi-inclusive bbh production rates that come 
from the uncertainties in the PDFs using the CTEQ PDF sets. First, the central 
value cross section (Tq is calculated using the global minimum PDF (i.e. CTEQ6M). 
The calculation of the cross section is then performed with the additional 40 sets 
of PDFs to produce 40 different predictions, cxj. For each of these, the deviation 
from the central value is calculated to be Aaf = |crj — ao\ when cTj^cxo. Finally, to 
obtain the uncertainties due to the PDFs the deviations are summed quadratically 

as Ao"^ = \f^~Aof^ and the cross section including the theoretical uncertainties 
arising from the PDFs is quoted as o^oHact-- 

In Figs. 7 and 8 we plot the total SM NLO cross section for bg bh at the 
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Figure 7: Comparison between theoretical uncertainties due to scale dependence and 
uncertainties arising from the PDFs at the Tevatron for semi-inclusive bh produc- 
tion in the Standard Model. In the bottom plot, both uncertainty bands have been 
normalized to the central value of the total cross section ctq- 



Tevatron and LHC respectively. Here, we compare the uncertainties from residual 
scale dependence and the PDFs both for the total cross section (left) and the total 
cross section normalized to the central value calculated with CTEQ6M (right). 

From Fig. 8 one can see that, at the LHC, the theoretical uncertainty is dominated 
by the residual scale dependence. Due to the large center of mass (c.o.m.) energy of 
the LHC, the gluons and bottom quarks in the initial state have small momentum 
fraction (x) values and, hence, small PDF uncertainties typically in the 5-10% range. 

In contrast, due to the smaller c.o.m. energy, the PDF uncertainties at the Teva- 
tron (Fig. 7) are comparable and even larger than the uncertainties due to residual 
scale dependence over the full Higgs mass range. The smaller c.o.m. energy results 
in higher-x gluons and bottom quarks in the initial state which corresponds to large 
PDF uncertainties in the 10-30% range. 
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Figure 8: Comparison between theoretical uncertainties due to scale dependence and 
uncertainties arising from the PDFs at the LHC for semi-inclusive bh production in the 
Standard Model. In the bottom plot, both uncertainty bands have been normalized 
to the central value of the total cross section ctq. 
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